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We compute the ground-state fidelity and various correlations to gauge the competition between 
different orders in two-dimensional t-J-type models. Using exact numerical diagonalization tech- 
niques, these quantities are examined for (i) the plain t-J and t-t'-J models, (ii) for the t-J model 
perturbed by infinite-range d-wave or extended-s-wave superconductivity inducing terms, and (iii) 
the t-J model, plain and with a d-wave perturbation, in the presence of non-magnetic quenched 
disorder. Various properties at low hole doping are contrasted with those at low electron filling. In 
the clean case, our results are consistent with previous work that concluded that the plain t-J model 
supports d-wave superconductivity. As a consequence of the strong correlations present in the low 
hole doping regime, we find that the magnitude of the d-wave condensate occupation is small even 
in the presence of large d-wave superconductivity inducing terms. In the dirty case, we show the 
robustness of the ground state in the strongly correlated regime against disorder. 
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I. INTRODUCTION 

Understanding the mechanism of high-temperature su- 
perconductivity has remained a subject of much inter- 
est since its experimental discovery in the cuprates in 
1986>i More recently, this subject has received renewed 
attention following the emergence of the first iron based 
(pnictide) high-temperature superconductor.— It is gen- 
erally believed that high-temperature superconductivity 
has its roots in the interplay of strong correlations and 
reduced dimensionalityi^i^i^ However, a full theoretical 
understanding of this phenomenon has proven challeng- 
ing, and consensus regarding its microscopic origin has 
not yet been reached i^i^ii^ii^ 

A further complication arises from experimental find- 
ings that the doped cuprates are highly inhomogeneousjiS 
This feature has been the subject of numerous recent ex- 
perimental studies using local probes such as scanning 
tunneling spectroscopy (STS) J^ii^ii^iiiii^ii^ Theoretical 
studies of correlations and disorder in superconducting 
lattice models have, in general, either focused on d-wave 
BCS phenomenology in the presence of impuritie a^'^'^^i^^ 
or on microscopic disordered t-J and Hubbard type mod- 
els, sometimes with the addition of short-ranged terms 
that favor superconductivity. One basic result of BCS 
phenomenology is that non-magnetic impurities suppress 
superconductivity more strongly in nodal systems, i.e., 
d-wave superconductors, than in conventional s-wave 
superconductorsJiii^iiSii^ This raises the question why 
superconductivity in the high-temperature cuprates ap- 
pears to be rather resistant to impurity disorder, al- 
though they have a d-wave superconducting order pa- 
rameter. A large part of the answer presumably involves 
the short coherence length, which is of the order of a 



few lattice constants in these systems, as opposed to the 
enormous values attained in conventional superconduc- 
tors. Another problem within the standard BCS phe- 
nomenology appears to be that the only way correlations 
enter is through superconducting pairing channels, ne- 
glecting potentially important effects due to the pres- 
ence of fluctuations towards other competing instabili- 
ties. The robustness of high-temperature superconduc- 
tivity and its density of states against disorder has been 
recently studied in the framework of Hubbard and t-J 
models using Bogoliubov-deGenne a^^'^^ and Gutzwiller 
mean-field theories;^ and exact diagonalizatioui^ 

In a spirit similar to previous studies^i^ in this work 
we use numerical diagonalization of finite clusters to ex- 
amine the effects of doping a strongly correlated Mott 
insulator within the t-J model. The t-J model can be 
justified microscopically either by a large U jt expansion 
of the one-band Hubbard model^ii^ or by a reduction of 
the three-band copper oxide model to an effective single- 
band modeli^i^^ The latter approach provides greater 
freedom for the allowed parameter ratio of J l\t\. Fol- 
lowing Ref. [2^, we also consider the t-J model with the 
addition of an infinite-range superconducting term. This 
term is tunable, and structured to induce either d-wave 
or extended-s-wavc superconductivity. Furthermore, we 
analyze the effects of quenched disorder in the ground 
state of these systems. 

In this work, we focus on three key observables that 
provide unique insights into the properties of the t-J 
model. The first of these observables is the ground- 
state fidelity metric g, defined below in Eq. PH)) . This 
quantity is related to the rate of change of the over- 
lap between the ground states of two Hamiltonians in- 
duced by a small change of a control parameter. The 
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ground-state fidelity, originally studied in the context 
of quantum information theory, has been shown to be 
a sensitive indicator of changes in the ground state of 
many-body systems, as they occur in quantum phase 
transitions,^i^i2^i^i^i^i^i2^i^ The other two observ- 
ables of interest we will study are the d-wave and 
extended- s- wave superconductivity condensate occupa- 
tions. They will be defined carefully in the next section. 

A few words motivating this project and the tools used 
are appropriate at this point. By changing control pa- 
rameters in a quantum many-body system, one may en- 
counter first-order or continuous quantum phase transi- 
tion. In finite systems, first-order transitions are easy 
to monitor, using the ground-state energy and density 
matrices. We have studied pair-pair (with d-wave and 
extended- s- wave symmetries) and density-dcnsity/spin- 
spin correlations. The former indicate superconducting 
instabilities and the latter indicate charge or spin order- 
ings. Continuous transitions are more subtle. Due to the 
Wigner-von Neumann non crossing rule, different states 
of the same symmetry approach each other near a tran- 
sition but do not cross; thereby leading to energy gaps. 
These energy gaps ultimately close for very large sys- 
tems at quantum phase transitions. In recent work, the 
fidelity [Eq. ([7])], and especially the associated fidelity 
metric [Eq. (fTUl) ]. have been shown to be sensitive indica- 
tors of continuous transitionS) ^^i^''i^^i^^i'^°i^^i'^^i'^'^ and we 
will see below that these also track level crossings quite 
well, as in general g exhibits a jump in those cases. 

The paper is organized as follows. First, the t-J model 
and the observables of interest are described in Sec. [TTl 
Sec. mil describes our study of the plain t-J and t-t'-J 
models when tuning the ratio J/t. In Sec. lIVl we examine 
the t-J model with the addition of d-wave and extended- 
s-wave superconductivity inducing terms. The effects of 
quenched disorder in the fidelity and superconductivity 
order parameters are analyzed in Sec. [V] In Sec. IVIl we 
conclude with a summary of our findings. 

II. PRELIMINARIES 
A. Model Hamiltonians 

The Hamiltonian for the plain t-J model can be written 

as 

+JJ2p(^^■^J-ln,n^P, (1) 

where c|s a-nd c^^ are the creation and annihilation op- 
erators for an electron with spin s =1,1 on a site i, 
~ Es ^Is'^is is the density operator, P is a projec- 
tion operator to ensure that there are no doubly occupied 
sites, i.e., it is assumed that the local Coulomb repulsion 



is very large such that two electrons (with antiparallel 
spin) cannot be on the same lattice site, and 

^^'^\'Y^c\^^ss'C,s' (2) 

ss' 

is the local spin operator [a are the spin-1/2 Pauli matri- 
ces). The sums in Eq. ([T]) run over nearest-neighbor 
sites. 

It has been discussed is several works that longer- 
range hoppings are needed to reproduce the Fermi sur- 
face and electron-hole asymmetries for different cuprate 
superconductorsj '^^'^^i'^''''^^i'^^'^°i"^^ Here we will briefly dis- 
cuss the effect of the next-nearest-ncighbor hopping 

Ur^-t' i^(4s. + H.c.)p, (3) 

where now the sum runs over next-nearest- 

neighbor sites. 

Within a number of mean-field theories r^i^i^ the 
phase diagram of the doped t-J model exhibits d-wave 
superconductivity. However, it is not clear how quan- 
tum fluctuations, which are expected to be significant in 
two dimensions, would affect such a state. As suggested 
by numerical diagonalizatio n^^i^^'^^ and quantum Monte- 
Carlo studies r^i^ competing states with other broken 
symmetries are expected to occur in the proximity of the 
superconducting phase. In particular, there has been a 
debate within the literature whether at low hole dop- 
ing, phase separation in the two-dimensional case oc- 
curs already at infinitesimal J/t > or only beyond 
a finite threshold value r^i^i^i^i^ preceded by d-wave 
superconductivity4^i^'^ It thus appears that the exact 
phase diagram of the t-J model still remains to be set- 
tled. Here, we address some of these issues using the 
fidelity as an indicator. Meanwhile, in order to precip- 
itate a superconducting ground state in the presence of 
strong correlations, we add attractive terms of the form 

n, = ~^Y.PVlv,P, (4) 

which favors a d-wave superconducting pairing, and 

ns = -^Y.Ps\s,P, (5) 

which favors an extended-s-wave superconducting pair- 
ing. 

In Eqs. (SI) and ®, A = (A^,»+x - A,,,+y), 
Si = (Ai_i+x -I- ^i,i+y) and K^j = Ci^Cji + Cj|Ci|. 
Ad and As denote the strengths of the d-wave and 
extended-s-wave superconductivity inducing terms, re- 
spectively, and L is the number of lattice sites. Equa- 
tions (jll) and ^ have infinite-range terms of the type 
Bardecn, Cooper, and Schrieffer considered in their re- 
duced Hamiltoniau)^ while building in the d-wave and 
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extended- s- wave symmetries in the superconductivity or- 
der parameter. Within mean-field theory, Eq. ^ with 
the addition of Eq. ^ leads to the same d-wavc ground 
state as found from the t-J model42ii^ 

Finally, we also consider the effects of quenched ran- 
dom disorder of the form 



random 



E 



(6) 



where the e^'s are taken randomly from a uniform distri- 
bution between [— r,r]. The full Hamiltonians given by 
Eqs. (HI)-® thus describe inhomogcneous strongly cor- 
related superconductors. In this study, we use numerical 
diagonalization of clusters with 18 and 20 sites and pe- 
riodic boundary conditions. The cluster geometries can 
be found in Ref. S The dimension of the largest Hilbert 
space diagonalized is of the order of 10^. 



B. Observables 

The first observable of interest is related to the fidelity 
F, which is defined as follows. Assume a general Hamil- 
tonian of the form 

W(A) = Tip + ATii, 

where Tii is taken to be the driving term. In the next 
sections, we will consider Tii to be cither the Heisenberg 
interaction term in Eq. (U), the superconducting terms 
given by Eqs. (j4]) and ([5|), or the disorder term in Eq. 
([6]). Let |4'o(A)) be the (normalized) ground state of 
7tI(A) and |^'o(A 4- (5A)) be the (normalized) ground state 
of 7i(A -I- 6\). The fidelity is then defined as the overlap 
between |*o(A)) and |*o(A + ^A)), i.e.. 



F{\5\) = |(*o(A)|«'o(A + a))|. 



(7) 



If the ground state is nondegenerate, and if 5X is suffi- 
ciently small, one can compute |5'o(A-|-(5A)) up to second 
order in perturbation theory. The only two terms of the 
(normalized) second order expansion that have a nonva- 
nishing overlap (n.o.) with |\I'o(A)) are 



|^o(A-H5A))„.o. = |*o(A)) X 



(8) 



|(vI/„(A)|7t:i|vI/o(A))|2 



where |5'q(A)) are the eigenstates of the Hamilto- 
nian with eigenenergies -Ea(A), i.e., ?T!(A)|*I'a(A)) = 
E^{X)\^^{\)). 

This means that up to the lowest order in 5\, one can 
write the fidelity in the form 

^(A,a)^l-^E^^^^^^^™^- (9) 



Since the sum on the rhs is in most cases an extensive 
quantity (see, e.g., Refs.[2l,[3l|, and [3^ for counterexam- 
ples see, e.g., Ref. 113), one can define the fidelity metric 

as 



5(A,<5A) 



lim q(\,5X) 



2 1-F(A,(5A) 
L 



\{^^{\)\H,\Mm' 



(10) 



In the following, we refer to lim5;j_»o g{X, 6X) as 5(A) or 
simply as g. From its definition, the fidelity metric g is 
dimensionless, positive and (in most cases) intensive, i.e., 
of 0(1). This is one of the main quantities that we will 
examine in the following sections. F{X,6X) will be com- 
puted using Lanczos diagonalization, choosing a value of 
i5A that is sufficiently small so that it does not affect the 
result of the ratio in Eq. pO)) . i.e., giving effectively the 
value in the lim5;)j^o 5(A, ^A). The above is of course true 
provided one does not encounter a level crossing. At a 
crossing, we compute g on either side of its jump by the 
above limiting process. 

The other two quantities of interest are the d-wave 
and extcnded-s-wave superconductivity condensate oc- 
cupations. Given the d-wavc pair density matrix 

Pfj^i^olPvjv.Pl'fo), (11) 
and the extcnded-s-wave pair density matrix 



(12) 



[EoiX) - E^iX)Y 



the d-wave (Af) and extended-s-wave (Af) condensate 
occupations are defined as the largest eigenvalues of Pfj 
and Pij^ The corresponding eigenvectors of the den- 
sity matrices are known as the "natural orbitals," and 
those with the largest eigenvalues are referred to as the 
"lowest natural orbitalsi^" If a condensate of pairs with 
a particular symmetry occurs in the system, the corre- 
sponding condensate occupation will scale linearly with 
the total number of fermions, as the system size L is 
increased while keeping the density constant 1^ This in 
turn is equivalent to stating that Pf^ and P^j exhibit off- 
diagonal long-range orderi^ Condensation also implies 
that all other eigenvalues are A^'* ~ 0(1) ,S An advan- 
tage of using these definitions is that they are valid in- 
dependently of whether the system is translationally in- 
variant or not, i.e., they work the same in clean systems 
and in the presence of disorder. In the particular case 
of translationally invariant systems, the eigenvalues Pf^ 
and P^j are occupations in momentum space. 

Since we will be dealing here with systems with differ- 
ent densities and finite sizes, in many cases we find it use- 
ful to monitor the ratios = Af/Aj and R" = Af/A| 
between the largest eigenvalues {Af, Af) and the sec- 
ond largest eigenvalues (Ag, Aj) of the density matrices. 
These ratios were first introduced in our earlier work^ 
and here we briefly reiterate the motivation behind this 
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construction. If condensation occurs, i.e., symmetry is 
broken in the thermodynamic limit, these are equivalent 
to studying Af and Af because the next eigenvalue is 
small, i.e., Aj"' ~ 0(1). However, computing i?'' and W 
has the added benefit of eliminating uninteresting nor- 
malization effects related to the change in the particle 
density, etc. It also has some advantages when trying to 
understand the effects of changing a Hamiltonian param- 
eter for a system with a fixed size, where we find cases 
with Ai and R behaving differently. 



III. PLAIN t-J MODEL 

As a first step, in this section we study how the ob- 
servables of interest behave within the plain t-J model. 
We begin with the effect of the antifcrromagnetic Heisen- 
berg coupHng on the ground state of this system. Within 
many mean-field theories, finite values of J favor su- 
perconducting ground states close to half-filling.— i^i^ 
The energy scale is set by the hopping parameter t = 1. 
We further study the ground state of the t-J model for 
values of J between and 1. While this rather large 
range is not achievable within the large U/t expansion 
of the one-band Hubbard model (where J ~ t'^/U), it 
should rather be regarded as a "Gedanken range," in- 
tended for the purpose of studying the possible phases 
that emerge, including the widely discussed possibility of 
phase separation4ii^i^ 

We begin by computing the fidelity metric in Eq. PU)) . 
considering the Heisenberg coupling in Eq. ([1]) to be the 
tuning parameter in the Hamiltonian, so that J plays 
the role of A. For our computations, we find that taking 
5 J = 10~^ is small enough to provide results for g which 
are independent of 6 J. 

In Fig. [TJ we compare the dependence of the fidelity 
metric on J for low electron filling [Figs.[T]Ja) and[ljb)] 
and low hole filling [Fig. [Ijc)]. The contrast between the 
two cases is very clear. For low electron filling the fidelity 
changes very little between J = and 1, but we do ob- 
serve some level crossings, signaled by relatively small 
jumps in the fidelity metric. These have their origin in 
the peculiarities of particular cluster shapes, as exempli- 
fied by the fact that for the same particle number they 
are found in one of the clusters and not in the other (see, 
e.g., the case of ten particles), and therefore they are 
expected to disappear in the thermodynamic limit. 

For low hole doping, on the other hand, see Fig. [TJc), 
the fidelity metric exhibits a large response for small 
J < 0.05, suggestive of a continuous phase transition 
in that regime. The cluster sizes accessible to numerical 
diagonalization are too small to identify a critical point, 
if any, but we find that a qualitatively similar behavior 
occurs for two and four holes in both clusters with 18 and 
20 sites. 

Next consider the occupation of the d-wave and 
extcnded-s-wave lowest natural orbitals as J is tuned in 
the system. Results for these quantities are shown in Fig. 



0.02 



0.01 



0.01 

0.005 


100 



(a) 



L=18 



■ N=4 
N=6 ■ 
N=8 
N=l6 



(b) 



L=20 



10 



oo 1 



0.1 



0.01 





L=18, N=14 '\ 




L=18, N=16 ' 




L=20, N=16 1 


I 


L=20, N=18 : 



















0.2 



0.4 



J 



0.6 



0.8 



FIG. 1: (Color online) Fidelity metric of the t-J model as 
a function of the antiferromagnetic coupling J. (a) Cluster 
with 18 sites and fillings N = 4,6, 8, 10. (b) Cluster with 20 
sites and fillings A'' = 4, 6, 8, 10 [the same parameter sets as 
in (a)], (c) Clusters with 18 and 20 sites and two and four 
holes. Notice that for all our results in this paper we consider 
N = N^+Ni, where At = (c|c|), Ni = {c|cj>, and At = A|. 



[5] As observed here, changes in the fidelity metric in Fig. 
[T] are accompanied by changes in the lowest natural or- 
bital occupations in Fig. [51 Interestingly, for each given 



0.6 
J3.4 
" 0.2 


0.4 



< 1.5 



- (a) ' ■ ' L=/S'. 






'■ N=iq_ 




--^-^ — ' — \ — ' — \ — 1 — 1 — 1 — 


(b) ' ' ' L=2o' ' 


(e) L=20 


— -■-^^ 




— 1 1 1 1 1 1 1 1 1 

-(c) .... 

' r " ^' ■ 


=^ — \ — < — \ — < — \ — 1 — 1 — 1 — 
- (f) 


-i L=18. N=14 . 

L=18.N=16 

-j L=20, N=16 - 

y. L=20, N=18 



1.5 < 



0.4 0.6 
J 



0.4 0.6 
J 



FIG. 2: (Color online) Occupation of the lowest natural or- 
bital, of the d-wave (left panels) and extended-s-wave (right 
panels) density matrices, as a function of J. [(a) and (d)] 
Cluster with 18 sites and fillings A = 4, 6, 8, 10. [(b) and (e)] 
Cluster with 20 sites and fillings A = 4, 6, 8, 10. Notice that 
the parameter sets for (a)-(e) are all the same and given in 
(d). [(c) and (f)] Clusters with 18 and 20 sites and two and 
four holes. The parameter sets in (c) and (f) are also the 
same. 
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cluster shape and filling fraction, the behavior of Af{J) 
and Af(J) is qualitatively similar. We do find, however, 
a marked difference between the results obtained for Af 
and between low electron filling, and low hole doping. 
For low electron filling [Figs, [^a), [^b), Hd), andllfe)], 
the lowest natural orbital occupation only changes ap- 
preciably when a level crossing is observed in the fidelity 
(Fig. [T|). Also, in this case, we find Af > A^^, suggest- 
ing dominant s-wave order parameter that is consistent 
with the results in Refs. [4^ and l47l. In contrast, for low 
hole doping [Figs.[2l[c) and[2tf)] with two and four holes 
in both cluster geometries, the occupation of the lowest 
natural orbitals in all cases exhibit a large increase for 
small values of J (J < 0.05), where a sizable response 
was observed in g in Fig. [T] 
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FIG. 3: (Color online) Ratio between the two lowest natu- 
ral orbitals, of the c?-wave (left panels) and extended- s-wave 
(right panels) density matrices, as a function of J. [(a) and 

(d) ] Cluster with 18 sites and fillings iV = 4, 6, 8, 10. [(b) and 

(e) ] Cluster with 20 sites and filhngs iV = 4, 6, 8, 10. Notice 
that the parameter sets for (a)-(e) are all the same and given 
in (a), [(c) and (f)] Clusters with 18 and 20 sites and two 
and four holes. The parameter sets in (c) and (f) are also the 
same. 



holes in 18 and 20 lattice sites, whereas R" decreases for 
the same set of parameters. This is a strong indication 
that for larger systems sizes the model may select, if any, 
a d-wave superconducting ground state. To further ex- 
plore this possibility, in the next section we perturb the 
t-J model Hamiltonian with d-wave and extcnded-s-wave 
superconductivity inducing terms and study its response. 

We have also studied other correlation functions, such 
as the spin-spin and density-density correlations. The 
results for the ratio between the two lowest natural or- 
bitals of those two matrices are presented in Fig. 0] for 
two and four holes in 18 and 20 sites. As expected, 
spin-spin correlations (/j«'*-co" -| exhibit a large enhance- 
ment corresponding to the onset of antiferromagnetic or- 
der. Jiss-corr j-gfj^^Qg dramatically as the hole doping 
increases, in contrast to the behavior of i?** in Fig. [3ljc), 
which for large J increases as the doping increases (in the 
underdoped regime). Interestingly, stripes are expected 
to be one of the competing orders at low doping, and they 
should be reflected in the density-density structure fac- 
tor (in Rdd-corry Fig.g^b), One can see that Rdd-corr 
actually decreases with increasing J, i.e., no signature of 
charge-density order can be seen in our clusters. 
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FIG. 4: (Color online) Ratio between the two lowest natural 
orbitals as a function of J in clusters with two and four holes 
with 18 and 20 lattice sites for (a) the spin-spin correlation 
matrix and (b) the density-density correlation matrix. 



It is interesting to observe that the occupation of the 
lowest natural orbitals of both d-wave and extended-s- 
wave symmetries in general increase as a result of in- 
creasing J. Since we cannot perform finite size scaling 
by studying larger systems, this does not throw light on 
the competition between these two orderings. In order to 
resolve this issue, we have also studied the ratio between 
the lowest two natural orbitals for both superconducting 
symmetries. The results are presented in Fig. |31 They 
show a quite different behavior from the one seen for the 
lowest natural orbitals in Fig. [51 We see from Figs. EJa), 
ISlJb), [njd), and[3]^e), that for low electron fiUing the ra- 
tio R between the two lowest natural orbitals in general 
remains almost unchanged when J takes values between 
and 1, for both d-wavc symmetry and extended- s-wave 
symmetry. For low hole doping, Figs.[3]Jc) and[3]Jf), the 
behavior is different. R'^ increases for both two and four 



A. Next-nearest-neighbor (f') term 

As mentioned in Scc.[Tll several authors have discussed 
the importance of introducing long-range hopping in mi- 
croscopic Hubbard and t-J-like models in order to cap- 
ture many of the features observed in high-temperature 
superconductorS i'^^'^^i'^^'^^i"^^'^^ Here, we briefly discuss 
how a finite next-nearest-neighbor hopping affects the 
results discussed above for the plain t-J model. 

In Fig. [5l we show the fidelity metric for exactly the 
same clusters and fillings as in Fig. [TJ The results in the 
presence of a finite and small t' are qualitatively simi- 
lar to those of the plain t-J model [Fig. [1] in both the 
low-electron- filling and low-hole-doping regimes, and for 
negative (left column) and positive (right column) val- 
ues of t! . It is worth noticing, however, that in the low 
doping regime the presence of a finite < [Fig. EJc)] 
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moves the maximum response of the fidelity metric to- 
wards larger values of J, while no such effect is seen for 
t' > [Fig. EJf)]. The former displacement for i' < is 
more pronounced for four holes than for two holes in both 
clusters with 18 and 20 sites. This is consistent with pre- 
vious works which have pointed out that in hole-doped 
systems t' competes with J as it suppresses both antifcr- 
romagnetic corrclations^SiSiffi and the superconducting 
Te )^°'^^ and that the effect of t' becomes stronger with 
doping (in the low- hole-doping regime). A positive t' , on 
the other hand, enhances antiferromagnetic correlations, 
which is consistent with the results in [Fig. [^f)] where 
the maximum response of the fidelity is in all cases seen 
at J = 0. 
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FIG. 5: (Color online) Fidelity metric of the t-t'-J model 
as a function of the antiferromagnetic coupling J for (a)-(c) 
t' = -0.3 and (d)-(f) t' = 0.3. [(a) and (d)] Cluster with 18 
sites and fillings iV = 4, 6, 8, 10. [(b) and (e)] Cluster with 20 
sites and fillings A'^ = 4,6,8, 10 [the legends in (a), (b), and 
(d) are the same as in (e)]. [(c) and (f)] Clusters with 18 and 
20 sites and two and four holes (the legends are the same in 
both figures). 

The results for the lowest natural orbitals and the ra- 
tio between the lowest the second lowest natural orbitals 
are also qualitatively similar to those of the plain t-J 
model, and will not be presented here. In the remaining 
of the paper, we will set t' = and study the effects of 
superconductivity enhancing terms and disorder in the 
t-J model. 



IV. SUPERCONDUCTIVITY INDUCING 
TERMS 

A. d-wave term 

Let us first consider a total Hamiltonian that is the sum 
of Eq. ID) and the d-wave-superconductivity inducing 
term in Eq. ([5]), and study the ground state of this model 



as a function of increasing the parameter A^. In the fol- 
lowing, we fix the Heisenberg coupling to be J = 0.3, 
which is a value commonly used in the t-J model litera- 
ture. From the analysis in the previous section, we know 
that, at least for the finite clusters considered here, no 
further qualitative changes occur in the observables of 
interest for larger values of J. 

Recall from previous work)^ that the added d-wave 
term Eq. (j4]), being of infinite range, must certainly pre- 
cipitate superconductivity in the d-wave channel. This 
is because mean-field theory becomes exact in the ther- 
modynamic limit for an infinite-range model of this type. 
(The same argument, of course, also works in the pres- 
ence of an extendcd-s-wave term.) While this argument 
is true for very large systems, for finite systems one may 
need a finite ~ 0(1) to achieve superconductivityi^ 
Therefore, we expect that in some cases the fidelity met- 
ric should show an enhancement as a function of A^ at 
some characteristic value AJ. A small value of AjJ, con- 
sistent with X*i ~ 0, may be taken as an indicator of the 
incipient order of the A^ = model (the plain t-J model). 

In Fig. [51 we show the fidelity metric as a function of 
the driving parameter A^. In our calculations, we have 
taken 6Xd = 10~^, which is sufficiently small to ensure 
results consistent with the limit 6Xd ^ 0. 
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FIG. 6: (Color online) Fidelity metric of the t-J model (J = 
0.3) after the addition of a d- wave-superconductivity inducing 
term, (a) Cluster with 20 sites and fillings iV = 4, 6, 8, 10. (b) 
Clusters with 18 and 20 sites and two and four holes. The 
inset in (a) depicts the results for the cluster with 18 sites 
and fillings = 4, 6, 8, 10. 

Results for low electron fillings are shown in Fig. [D^a) 
and its inset. In all cases one can see that there is al- 
most no response in g when A^ is small and that a strong 
response occurs for Aj 1, indicative of a phase transi- 
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tion for a finite value of \d- These results are consistent 
with the behavior of the ratios and which are 
depicted in Figs. (Ha), Wih), [^d), and W[c). For most 
low fillings, both ratios change very little for small val- 
ues of Xd. Around ~ 1, they either jump abruptly 
(cluster with L = 18) or increase rapidly (cluster with 
L = 20). Notice that for large Ad there is almost one 
order of magnitude difference between the ratios seen for 
the occupation of the d-wave related natural orbitals and 
the extended-s-wave related orbitals. This is expected 
since the driving term has d-wave symmetry and hence 
d-wave superconductivity should be stabilized for large 
values of A^. 

The results for low hole doping (two and four holes) are 
in contrast with those of low electron filling. Figure El^b) 
shows that in the former case g exhibits a large response 
for very small values of A^. The behavior of g in this case 
is consistent with a phase transition at Xd ~ 0. The situ- 
ation is similar to that of g in the one-dimensional Hub- 
bard model as one tunes the onsite repulsion parameter 
U where the Mott phase transition occurs ai U = 0. 
In addition, as shown in Fig.[7]Jc) and[7Kf), the response 
of g for small values of A^ is accompanied by a continuous 
increase of R"^ and a continuous decrease of i?* for small 
values of A^. 




FIG. 7: (Color online) Ratio between the two lowest natu- 
ral orbitals, of the d-wave (left panels) and extended-s-wave 
(right panels) density matrices, as a function of A^. [(a) and 

(d) ] Cluster with 18 sites and fillings iV = 4, 6, 8, 10. [(b) and 

(e) ] Cluster with 20 sites and fillings iV = 4, 6, 8, 10. Notice 
that the parameter sets for (a)-(e) are all the same and given 
in (d). [(c) and (f)] Clusters with 18 and 20 sites and two 
and four holes. The parameter sets in (c) and (f) are also the 
same. 

Comparing the results in this subsection with Sec. IIIIl 
we find support for the view that in the thermodynamic 
limit, the plain t-J model is superconducting (with d- 
wave symmetry, for finite values of J), without the need 
of introducing A^. Finite values of A^ certainly enhance 
the superconducting features of the t-J model in finite 
clusters but may not be needed for larger system sizes. 



Earlier evidence in this direction comes from high tem- 
perature expansion studies^ and exact diagonalization 
studies of the plain t-J modelM^^^^ 

We should stress that the magnitudes of the ratios R"^ 
in Fig. Wic), reveals a very important characteristic of 
the d-wave superconducting state generated in a lightly 
doped Mott insulator. The d-wavc condensate occupa- 
tion is low if one compares it with the one generated, 
for low electron filling [Figs. [Tl^a) and[7^b)], under the 
influence of a suitable value of A^- This low condensate 
occupation is a direct consequence of the strong correla- 
tions present in the system and is similar to the behavior 
of the condensate occupation in liquid helium, which is 
known to be strongly depleted due to the effects of strong 
interactions. 



B. s-wave term 

In order to further study the suggestions in Sec. IIV Al 
we have also studied the effect of adding an extended- 
s-wave-superconductivity inducing term [Eq. ^] to the 
plain t-J model. In Fig. [HI we explore the behavior of 
the fidelity metric for such a model. 
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FIG. 8: (Color online) Fidelity metric of the t-J model 
(J = 0.3) after the addition of an extended-s-wave- 
superconductivity inducing term, (a) Cluster with 20 sites 
and fillings iV = 4, 6, 8, 10. (b) Clusters with 18 and 20 sites 
and two and four holes. The inset in (a) depicts the results 
for the cluster with 18 sites and fillings A'^ = 4, 6, 8, 10. 

Figure [5] shows that in this case small values of the 
control parameter As < 1 induce almost no response in 
the fidelity metric at any filling. Instead, one needs a 
sizable magnitude of Ag, at both low electron filling and 
low hole doping, to trigger a large response in g, which 
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would signal a phase transition to an extended-s-wave 
superconducting state. The calculated ratio between the 
two lowest natural orbitals, depicted in Fig. [9l is also 
consistent with this observation. Only large values of 
As enhance the extended-s-wave condensate occupation 
for all fillings [Figs. [9l^d)-[9l^f)], particularly for low hole 
doping [Fig.[ni^f)]. Interestingly, in analogy to the d-wave 
case for W^, the enhancement of the ratio is almost 
an order of magnitude larger at low electron fillings than 
at low hole doping, evidencing the strong depletion of 
the condensate that occurs in the latter case due to the 
presence of strong correlations. 




FIG. 9: (Color online) Ratio between the two lowest natu- 
ral orbitals, of the d-wave (left panels) and extended-s-wave 
(right panels) density matrices, as a function of As. [(a) and 

(d) ] Cluster with 18 sites and fillings iV = 4, 6, 8, 10. [(b) and 

(e) ] Cluster with 20 sites and fillings A'' = 4, 6, 8, 10. Notice 
that the parameter sets for (a)-(e) are all the same and given 
in (d). [(c) and (f)] Clusters with 18 and 20 sites and two 
and four holes. The parameter sets in (c) and (f) are also the 
same. 

Overall, these findings are consistent with previous ex- 
act diagonalization studies of the t-J model^i^i^ that 
have provided numerical evidence for dominant d-wave 
pairing at low hole doping, absent in the opposite limit 
of low electron filling. Our results also show that strong 
correlations become more important as one approaches 
half-filling, where no matter which superconductivity en- 
hancing term one introduces, the resulting condensate 
occupation is always small. 

V. QUENCHED DISORDER 

In this section, we address the question of the effect 
of disorder in the models and quantities studied in the 
previous sections for the clean case. 

For the plain t-J model, in Sec. Illli we have seen 
that the lowest natural orbital and the ratio between the 
lowest and second lowest natural orbitals are in general 
small. For that model, we focus here on the response 



of the fidelity metric to adding disorder to the system. 
Results for this quantity, averaged over different disorder 
realizations, are presented in Fig. [TOl For weak inter- 
actions, one would expect small values of F to produce 
localization and dramatically change the nature of the 
ground state in the clean case. In contrast, for strong 
correlations, this figure illustrates that, except for the 
very low electron-filling case of four particles, one needs 
a very large disorder strength (F > 1) in order to trig- 
ger a large response in g. Different disorder realizations 
produce a large response in g for different values of F, 
hence, the various peaks that can be seen in many of 
the plots (notice that the results for g are presented in 
a semi- log scale). However, the common feature for all 
fillings (with the exception of = 4) is that large re- 
sponses only occur when F > 1. This indicates that the 
nature of the ground state of the system in these cases 
is robust against nonmagnetic disorder and exemplifies 
the importance of strong correlations in the system. The 
exception is the case of very low fillings [e.g., = 4 for 
L = 18 and L = 20 in Fig.finia) andfTn^b), respectively] 
in the t-J model, where the Gutzwillcr projection plays 
a subdominant role. 
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FIG. 10: (Color online) Fidelity metric of the t-J model (J = 
0.3) as a function of the disorder strength F. (a) Cluster 
with 18 sites and fiUings iV = 4, 6, 8, 10. (b) Cluster with 20 
sites and fillings TV = 4, 6, 8, 10 [the same parameter sets as in 
(a)], (c) Clusters with 18 and 20 sites and two and four holes. 
In all cases the fidelity metric was computed averaging over 
20 different disorder realizations, except for the cluster with 
L — 20 and fillings N = 10 and N — 16 where the average 
was performed over ten disorder realizations. 

In the clean case, the addition of a d-wave supercon- 
ductivity inducing term was used to precipitate d-wave 
superconductivity in the t-J model and was shown to pro- 
duce a large enhancement of the ratio R'^ for low electron 
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fillings and a smaller enhancement for low hole doping. 
Here, we study the effect of disorder in both the fidelity 
metric g and the ratio i?'^ when J — 0.3 and = 1.25. 

In Fig. [Til we have plotted the evolution of g (left 
panel) and R'^ (right panel) with increasing disorder. The 
results depicted in this figure were obtained averaging 
over the same number of disorder realizations as in Fig. 
[TOl By comparing g in Figs. [TOl andfTTl one can see that 
the effect of adding a finite value of \d ~ 1.25 is more 
pronounced at low fillings, where large responses in g 
move toward larger values of F (F > 2, with the exception 
of = 6 in the i = 18 cluster that still has a peak 
around F = 1). These large responses are expected to 
indicate of the destruction of the superconducting ground 
state generated in the clean case. This can be better seen 
in Figs. [TIT d) and [TlTe). where the value of R'^ is much 
smaller than the one in the clean case whenever a large 
response is seen in g [Figs. [TIJa) and fTTT b)]. 

Interestingly, in the low hole-doping regime, the be- 
havior of the fidelity metric as a function of disorder is 
very similar when = 1.25 and in the clean case, an 
indication that in this regime the ground state of the t-J 
model is not dramatically affected by the presence of Ad. 
This marked contrast with the low electron filling is fur- 
ther supported if one realizes that the relative reduction 
of i?*^ is much more drastic at low electron filling than in 
the low hole-doping regime. Further studies using alter- 
native numerical approaches may be needed in to clarify 
the scaling of these effects with system size. 

Overall, our finding in the presence or absence of the 
superconductivity inducing term is that strong correla- 
tions generate a groimd state in the t-J model that is 
robust against disorder. 



VI. SUMMARY 

Within the t-J model, we have calculated the ground- 
state fidelity, the d-wavc and extended-s-wave condensate 
occupations, with and without the addition of a super- 
conductivity inducing term, and in the clean and disor- 
dered cases. 

In the clean case we find that: 

(i) the plain t-J model exhibits a distinctive signature 
in the fidelity metric at low hole doping and small val- 
ues of J that may be an indication that a continuous 
phase transition occurs in this regime. By studying the 
d-wave and extended-s-wave condensate occupations, we 
find that the former is favored. As expected, spin-spin 
(antiferromagnetic) correlations are enhanced by J in the 
low hole-doping regime, but they are drastically reduced 
by doping. Density-density correlations, on the other 
hand, do not exhibit any clear signature of order as J is 
increased. Adding a i' < term only moves the response 
of the fidelity metric towards larger values of J. This is 
indicative of a continuous phase transition for a value of 
J that is larger than for the plain t-J model. 

(ii) In order to better understand the nature of the 
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FIG. 11: (Color online) Fidelity metric (left column) and ratio 
between the two largest eigenvalues of the d-wave density ma- 
trix (right column) for the t-J model (J = 0.3) with the addi- 
tion of a d-wave superconductivity inducing term (A^ = 1.25) 
as a function of the disorder strength P. [(a) and (d)] Cluster 
with 18 sites and fillings iV = 4, 6, 8, 10. [(b) and (e)] Cluster 
with 20 sites and fillings A'^ = 4, 6, 8, 10 [the same parameter 
sets as in (a) and (d)]. [(c) and (f)] Clusters with 18 and 20 
sites and two and four holes. In all cases g and R"^ ware com- 
puted averaging over 20 different disorder realizations, except 
for the cluster with L = 20 and fillings N = IQ and iV = 16 
where the average was performed over ten disorder realiza- 
tions. 



phases in the plain t-J model, we added infinite range 
d-wave and extended-s-wave superconductivity inducing 
terms. For a finite value of J = 0.3, we find that at low- 
hole doping arbitrarily small d-wave driving terms induce 
a large response in the fidelity metric, which is consistent 
with the plain t-J model having a d-wavc superconduct- 
ing ground state. In contrast, almost no response was 
found for low electron filling, and for the extended-s-wave 
superconductivity driving term for all fillings. In these 
latter cases one needs a large value of the driving term in 
order to obtain a sizable response. Interestingly, at low 
hole doping, we always find the condensate occupations 
to be small, i.e., there is a strong depletion of the (rele- 
vant) condensate no matter the symmetry of the driving 
term. This is an indication that, for the plain t-J model, 
mean-field theories based on the d-wave order parameter 
may not be justified. 

In the dirty case, we have shown that very strong dis- 
order is required to produce large changes in the ground 
state of the plain t-J model whenever the electron fill- 
ing is not very low, i.e., whenever the system is strongly 
correlated. Adding a d-wave superconductivity inducing 
term was shown to leave the results for the low hole- 
doping regime almost unchanged with respect to the ones 
of the plain t-J model, while the results for the low 
electron-filling regime were modified to a larger extent. 
This is, once again, consistent with the hypothesis that 
the plain t-J model has a d-wave superconducting ground 
state for low hole doping. Future studies on the scaling 
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of the effects and observables discussed here, with sys- 
tem size, could shed further Hght on the nature of the 
ground state of the t-J model and its relation to high- 
temperature superconductivity. 
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Note added on proof. -With regard to the ratios (con- 
densate occupations) R'^ and i?" in the low hole doping 
regime, we would like to emphasize that in addition to 
being small (even in the presence of large superconduc- 



tivity enhancing terms in the Hamiltonian) they increase 
linearly with increasing doping. This is shown in Fig. 1121 
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FIG. 12: (Color online) Ratio between the two lowest natural 
orbitals, of the d-wave (a) and extended-s-wave (b) density 
matrices, as a function of hole doping {uh) for a large su- 
perconductivity inducing term with the same symmetry of 
the order parameter plotted. Results are presented for clus- 
ters with 18 and 20 sites. For low hole doping one can see a 
nearly linear dependence of the ratios and W with Uh- 
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